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THE PHYSICAL CONDITIONS FOR MASSIVE STAR FORMATION: 
CONTINUUM MAPS AND MODELING 
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ABSTRACT 

Fifty-one dense cores associated with water masers were mapped at 350 /zm. These cores are very 
luminous, 10'^ < Lboz/L© < 10^, indicative of the formation of massive stars. Dust continuum contour 
maps, radial intensity profiles, and photometry are presented for these sources. The submillimeter dust 
emission peak is, on average, nearly coincident with the water maser position. The spectral energy 
distributions and normalized radial profiles of dust continuum emission were modeled for 31 sources 
using a one-dimensional dust radiative transfer code, assuming a power law density distribution in the 
envelope, n = nf{r/rf)~P. The best fit density power law exponent, p, ranged from 0.75 to 2.5 with 
(p) — 1.8 ± 0.4, similar to the mean value found by Beuther et al. (2002) in a large sample of massive 
star forming regions. The mean value of p is also comparable to that found in regions forming only 
low mass stars, but {uf) is over two orders of magnitude greater for the massive cores. The mean p is 
incompatible with a logatropic sphere (p = 1), but other star formation models cannot be ruled out. 
Different mass estimates are compared and mean masses of gas and dust are reported within a half-power 
radius determined from the dust emission, (log(M(< rdec))) = 2.0 ± 0.6, and within a radius where the 
total density exceeds 10^ cm~^, (log(M(< r„))) = 2.5 ±0.6. Evolutionary indicators commonly used for 
low mass star formation, such as T^oi and Lboi/ Lsmm, may have some utility for regions forming massive 
stars. Additionally, for comparison with extragalactic star formation studies, the luminosity to dust 
mass ratio is calculated for these sources, {Lboi/Mo) — 1.4 x lO'' Lq/Mq, with a method most parallel 
to that used in studies of distant galaxies. This ratio is similar to that seen in high redshift starburst 
galaxies. 

Subject headings: ISM: dust — stars: formation, high-mass — submillimeter 



1. INTRODUCTION 

The study of regions forming massive stars is essen- 
tial to our understanding of how stars are born. Since 
most stars form in clusters associated with high mass stars 
(e.g.. Carpenter 2000), many recent studies have focused 
on better understanding the physical conditions in these 
regions (e.g., van der Tak et al. 2000, Sridharan et al. 2002, 
Beuther et al. 2002, HatcheU et al. 2000, Walsh et al. 2001, 
Osorio, Lizano, & D'Alessio 1999, Garay & Lizano 1999). 
The density distribution in the envelopes of regions form- 
ing massive stars is an important observational constraint 
for star formation models. The density distribution is usu- 
ally a power law (n cx r~^). For example, McLaughlin & 
Pudriz's (1997) logatropic sphere model predicts a shallow 
power law with p = 1 in the static envelope, whereas Shu's 
(1977) inside-out collapse model for isolated star formation 
has an initial density distribution with p = 2. 

A power law density distribution has been fitted to 
observations of low mass star forming regions. Shirley 
et al. (2002a) and Young et al. (2002) find a combined 
(p) = 1.6 ± 0.4 for Class (Andre, Ward-Thompson, & 
Barsony 1993) and Class I sources (Lada & Wilking 1984, 
Lada 1987, Myers & Ladd 1993, Chen et al. 1995). All "er- 
rors" on mean values in this paper refer to the standard 
deviation of the distribution of values about the mean. 
Shirley et al. (2002a) and Young et al. (2002) also report 
that aspherical cores have shallower power laws. If the as- 
pherical cores are left out of the average for the low mass 



cores, (p) goes to 1.8. While low mass cores are well- 
studied, it is only recently that the density structure of 
high mass cores has been investigated for large samples. 
For example, van der Tak et al. (2000) found a shallow 
density structure, (p) = 1 to 1.5, for a sample of 14 re- 
gions forming massive stars, while Beuther et al. (2002) 
reported (p) = 1.6 ± 0.5 for a larger sample of 69 massive 
star forming regions. 

Studies of massive star forming cores also have impor- 
tant implications for understanding extragalactic star for- 
mation, including starburst galaxies. The far-infrared lu- 
minosity to dust mass ratio, i/M, is a tool often used in 
extragalactic studies to characterize star formation since 
it is proportional to the star formation rate per unit mass 
(Kennicutt 1998). To learn if starburst galaxies are form- 
ing stars by mechanisms similar to those in the Milky 
Way, but on grander scales, it is important to investigate 
the star formation efficiency and L/M for dense gas in 
more accessible Galactic star formation regions to provide 
a point of comparison between these modes of star forma- 
tion. 

1.1. The Sample 

The objects in this study were selected from the sample 
of Plume et al. (1992, 1997) of massive star forming cores 
associated with water masers. Table 1 lists the sources 
and their observed properties. Water masers are asso- 
ciated with regions of very dense gas {n > 10^" cm~^; 
Elitzur et al. 1989). Each of the cores had been mapped 
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in the CS J = 5 -> 4 transition (Shirley ct al. 2002b) 
and detected in the CS J = 7 — > 6 transition (T4 > 1.0 K; 
Plume et al. 1992). The critical density of CS J = 5 ^ 4 is 
ric = 8.9 X 10^ cm~^; however, a density, Ueff, of 2.2 x 10^ 
cm^'' will produce an observable line of 1 K for a gas tem- 
perature of 10 K (Evans 1999). For a gas temperature of 
100 K, which may better describe massive cores, the effec- 
tive critical density is even lower, rie// = 6.0 x 10"* cm~^ 
(Evans 1999). Consequently, models constrained by mul- 
tiple transitions arc needed to determine density. Plume 
et al. (1997) reported (log(n)) = 5.9 from LVG models 
of multiple CS transitions for the regions from which our 
sample was taken. Therefore, these objects were known to 
contain a significant amount of dense gas; however, their 
mass and density structures were not well known. Many 
(43%) of the regions in our study were also known to be 
associated with UCH II regions. 

The sample covers a large range of distances, from Ori- 
IRC2 at a distance of 450 pc to G12.21-0.10 at 13.7 kpc. 
The distances were found in the literature (see Table 1), 
and spectrophotometric distances were used when avail- 
able. The mean distance in the sample is 3.9 kpc, while 
the median distance is 2.8 kpc. Figure 1 shows the distri- 
bution of distances in the sample. 

Our sample includes nine objects described in van der 
Tak et al. (2000). The van der Tak et al. (2000) sources 
were selected to be luminous, visible from the North- 
ern hemisphere, and, in most cases, bright in the mid- 
infrared. For comparison, the sample of high mass proto- 
stellar objects studied by Beuther et al. (2002) were se- 
lected from objects north of —20° declination detected in 
CS J = 2 ^ 1 (Bronfman et al. 1996), with far- infrared 
colors characteristic of UCH II regions. Their sources 
are also bright at far-infrared wavelengths. However, the 
sources in the Beuther et al. (2002; also Sridharan et al. 
2002) sample were also chosen to have low radio contin- 
uum flux (< 25 mjy; Sridharan et al. 2002) to ensure that 
their sources were isolated and not typical UCH II regions. 
Surprisingly, our sample has only two objects in common 
with that of Beuther et al. (S231 and G12.89-F0.49). 

2. OBSERVATIONS 

2.1. Observation Techniques 

Fifty-one regions forming massive stars were observed 
with SHARC (the Submillimeter High Angular Resolu- 
tion Camera), described by Hunter, Benford, & Serabyn 
(1996), during 5 nights in 1997 (December 21 and 22) and 
1998 (July 15, 23, and 25) on the 10.4 m Caltech Submil- 
limeter Telescope. SHARC is a one-dimensional bolometer 
array with a FWHM beam size, 6mb, of 14" (see Figure 2 
and §2.3). 

SHARC 's linear array consists of 24 detectors (Hunter 
et al. 1996); therefore, the telescope must be scanned in 
azimuth at constant elevation to map the source. Each 
350 fim map consists of approximately 11 scans extend- 
ing 240" in azimuth scanning at a rate of 4" per second. 
The individual scans are shifted by 4" in elevation to ex- 
tend the mapped region and to eliminate gaps in the map 
due to bad pixels (pixels 1, 5, 15, and 16 of the 24 pixels 
in the SHARC array). The scanning rate and elevation 
shifts were selected to be slightly smaller than the size of 
the pixel, 5" in the focal plane of the array, to obtain 



better sampling (Hunter et al. 1996). The secondary was 
chopped at 1.123 Hz with a chop throw of 90" to 100" 
in the azimuth direction. SHARC observations were con- 
ducted only during very dry conditions with Tcso < 0.06 
(see §2.3 below). 

2.2. Image Reduction 

All of the data were reduced and restored with the 
standard program CAMERA. The restoration algorithm 
is based on a technique described by Emerson, Klein, & 
Haslam (1979). The despiking routine was used on those 
maps that had pixels with spikes above 10 a, which we 
identified by visual inspection during reduction. The rou- 
tine replaces the spiked pixel with the average value of 
the adjacent pixels. In some cases the source was highly 
peaked, so a higlica' sigma was used to ensure the cen- 
tral pixel was not removed by the despiking routine. The 
night of 1998 July 15 was unusually windy. Five maps 
made on this night (W28A2, G12. 89+0.49, G12. 21-0.10, 
G24.49-0.04, W43S) showed signs of being affected by the 
wind and were corrected with linear destriping. The de- 
striping affected the maximum pixel value, in most cases 
decreasing it by less than 10%. 

After the data were reduced with CAMERA, gray scale 
images of the restored and combined bolometer maps 
were created. The Image Reduction and Analysis Facil- 
ity (IRAF) was used to find the value of the maximum 
pixel, the average background, and a for each map. The 
voltages in each map were multiplied by the extinction 
correction, 

Vcorr ~ ^obs^ ^ ^ (1) 

where sec(z) is the average airmass at the time of obser- 
vation. The determination of r is described below in §2.3. 

Figures 3 7 show the dust continuum contour maps. 
The contour levels are even multiples of a or 10% or 20% 
of the peak signal with an additional lowest contour at 3 a. 
The (0,0) position is the location of the water maser from 
Plume et al. (1992, 1997). The positions of known UCH II 
regions are indicated by plus signs on the contour maps. 
Three of the 22 marked UCH II regions (G12.21-0.10, 
G23.95-I-0.16, and W43S) were listed as UCH II regions by 
Wood & Churchwell (1989) but, in the same study, were 
reported to have diameters greater than 0.1 pc. There- 
fore, these sources may also be classified at compact H II 
regions. 

Normalized radial intensity profiles were created as in 
Shirley et al. (2000). The intensity was azimuthally aver- 
aged and normalized to the peak emission. The normalized 
intensities, 7(6) //(O), were plotted versus the impact pa- 
rameter, b = 9D, a line of sight offset from the center by 
an angle 6, for a source at distance D. The radial profiles 
are truncated at a radius, Vprof, when the signal fell to 1 
a or at 60", whichever is smaller. We do not considerv- 
data beyond a 60" radius where simulation in our models 
of the effects of chopping becomes problematic. Photom- 
etry was also provided by the radial profile program. The 
sky-subtracted fiuxes for 30" and 120" diameter apertures, 
9ap, are listed in Table 1. 

2.3. Calibration 
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The extinction coefficients at 350 /zm, T350, were de- 
termined using skydips from the CSO tipper (tcso) at 225 
GHz. A scaling between Tcso and T350 has been determined 
by comparing skydips between the 225 GHz and 350 /im 
tippers, 

r350 = (23.5 ± 0.2)tcso (2) 

(R. Chamberlin 2000, private communication). To check 

the relationship between T350 and Tcso, Uranus was ob- 
served as it set during two nights in 1998 July with ex- 
ceptionally stable sky opacity (ar/ir) < 0.06). The vari- 
ation of peak voltage on Uranus with airmass determined 
T350. The resulting r35o was consistent with Equation (2). 
Since wc were finable to observe a source as it set during 
each night, and since rsso may vary throughout the night, 
we determined the extinction correction, Equation (1), for 
each image by scaling from Tcso measurements using Equa- 
tion (2). 

Full maps of Uranus and secondary calibrator sources, 
NGC 2071IR and W3(0H), were used to determine the cal- 
ibration factors, C*, for each observing run. Sky subtrac- 
tions were made for each image by measuring the voltage 
(Vsfcj,) through multiple 20" apertures {9 sky) away from 
the source and averaging the measured sky voltage. 

Total voltages measured in apertures of diameter 6ap 
were then corrected for sky emission by 

(3) 

The calibration factors were calculated for images 
taken in the 1998 July run by measuring the flux at 350 /im 
in 30" and 120" diameter apertures of two maps of Uranus. 
The total flux of Uranus in 1998 July was 266.5 Jy. Cal- 
ibration factors for sources observed in 1998 July were, 
with statistical errors, = (9.6 ± 1.5) mJy and 
= (8.5 ± 2.0) mJy for 30" and 120" apertures, 
respectively. Because no planets were visible to serve as 
pointlike sources in 1997 December, maps of NGC 2071IR 
and W3(0H) were used for calibration. The assumed 
fluxes for NGC 2071IR and W3(0H) were 177 Jy beam"! 
and 498 Jy beam~^, respectively (Sandell 1994). Volt- 
ages measured in 30" and 120" apertures in three maps 
of each secondary calibrator were used to determine cal- 
ibration factors. The 1997 December calibration factors 
averaged over the six maps were, with statistical errors, 
C^o = (9.6±1.3) mJy V"! and C^^o = (4.7±0.7) mJy 
While some of the statistical uncertainties are smaller than 
20%, wc expect systematic errors to not allow fluxes to be 
determined at 350 /im to better than 20% (Hunter et al. 
2000). Furthermore, calibrations in 1997 December deter- 
mined from NGC 2071IR and W3(0H) are expected to be 
worse since neither calibrator is a point source (see Fig. 
3 for a map of W3(0H)). This effect is obvious in C^^o, 
which is nearly a factor of two lower when using the sec- 
ondary calibrators. C^" was not affected by the extended 
emission and is stable from 1997 December to 1998 July. 
For this reason, C^^° from 1998 July is used to calibrate 
all of the 120" aperture fluxes. 

Radial beam profiles were obtained from two maps of 
Uranus in 1998 July (Figure 2). The data were binned in 
5" bins. Broad sidelobes are seen beyond 15". At 40", the 



sidelobe power is as high as —17 dB. The general shape 
of the two profiles is consistent, but there are variations 
in the strength of the sidelobes. Fitting the beam profiles 
with a Gaussian, we find that the FWHM is at least 14" 
which is significantly different from the beam size reported 
in previous studies using SHARC at 350 /xm (10", van der 
Tak et al. 2000; 11", Hunter et al. 2000). The difference 
could be the result of poor focus or temporal changes in 
the beam. Young et al. (2002) report a 3" change in the 
FWHM of the James Clerk Maxwell Telescope beam at 
450 /xm in the course of a single night. The beam mea- 
surements were not frequent enough to characterize the 
changes of the CSO beam over the course of the observa- 
tions in this study; therefore, we adopt 9mb = 14". 

3. RESULTS 

Most of the cores appear slightly elongated in the con- 
tour maps (Figures 3 7). This is likely an effect of chop- 
ping during observations as the maps are usually elon- 
gated along the chop direction. The chop direction is in- 
dicated on the maps with an arrow. However, some of the 
sources have extended asymmetrical emission at 350 /im 
distinct from the chopping asymmetry with an intensity 
of a few cr (e.g., 0N2S, G40.50-h2.54). Some (14%) of the 
maps have double or multiple peaks, indicating the pres- 
ence of more than one luminosity source. These sources 
are noted in Table 1. A few have close embedded double 
peaks (e.g., G23.95+0.16, S235), while others show more 
spatially-distinct peaks. For example, the NGC 7538 re- 
gion shows 3 distinct peaks in Figure 7; these were mapped 
separately, but are plotted in their relative locations. 

The contom maps also show that the 350 /im dust peak 
is often nearly coincident with the water maser position 
at the center of the map. Thirty cores (59%) have the 
dust centroid within Omb/"^ of the maser position (see Ta- 
ble 1). The mean absolute distance of the dust centroid 
to the maser position is 8". About half of the 22 known 
UCH II regions in the sample are also close (< Omb/"^) 
to the maser, however, the mean absolute distance of the 
UCH II position from the water maser is 10". This separa- 
tion is similar to that of the UCH II region from the dust 
peak with mean absolute distance of 11". Samples based 
on water masers favor an earlier phase of star formation 
than UCH II region samples (Cesaroni et al. 1988, Shirley 
et al. 2002b). The high coincidence of the dust centroid 
and maser position in these regions implies that the dust 
emission may be primarily tracing the earlier stages as 
well. 

The FWHM size of each source, 6 dec, was determined by 
deconvolving the telescope main beam from the observed 
FWHM of the core by subtracting 9mb from the observed 
FWHM in quadrature. The observed FWHM was deter- 
mined from the radial profile of each source. The decon- 
volved half-power radius is defined to be r^ec = {D/2)9dec, 
where D is the distance to the source. For the entire sam- 
ple, (rdec) = 0.16 ± 0.10 pc, with a median of 0.14 pc. 
Figure 8 shows the distribution of rdec- 

Photometry from the literature is collected in Table 
2. The observed spectral energy distribution (SED) was 
used to calculate the bolometric luminosity, Lboi , for each 
source, using only data taken with beam sizes > 20", with 
the exception of the endpoints to the SED. Luminosities 
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range over three orders of magnitude, about 10"^ to 10^ Lq 
with {Lboi) = 2.5 X 10^ L© and a median Lf,oi = 5.1 x 10** 
Lq . If all the luminosity is attributed to a single star, the 
range of spectral types would be B3 to 04, with the me- 
dian being an 09 star. Considering multiple sources and 
accretion luminosity would lower these spectral types. 

Many recent studies have used power laws to fit directly 
the radial intensity profiles and infer the density structure 
of star forming regions (e.g., Shirley et al. 2000, Beuther et 
al. 2002). Wc did not use this technique; however, for com- 
parison purposes, we describe some general trends in the 
radial profiles of our sample. Generally, the radial profiles 
follow a power law from about 12" (approximately 9mb) 
to 40". The profiles arc flattened toward the interior, and 
beyond 40", they sometimes deviate from the power law. 
The flattening at small angles is likely due to beam ef- 
fects, but could also result from fragmentation of the core. 
The change in the slope beyond 40" is not consistently 
steeper or shallower. Some radial profiles are distorted by 
the presence of multiple peaks. 

4. MODELS 

One of the major motivations for this study was to 
learn what density distributions fit the data. A power 
law density distribution was assumed of the form n(r) = 
nf{r/rf)~P, for p in a range of 0.5 to 2.5. We chose 
Tf = 1000 AU for convenient comparison to other stud- 
ies, but this rj is well inside our beam. The values of n/ 
should be taken only as indicative of likely mean densities 
as substructure is very likely on those scales. The observed 
radial intensity profiles and SEDs were modeled using a 
modified version of the one-dimensional dust continuum 
radiative transfer code by Egan, Leung, and Spagna (1988) 
and an observation simulation code described by Evans et 
al. (2001). The radiative transfer code calculates the ra- 
dial temperature distribution, T£,{r) self-consistently for 
each input model of n{r). 

We also included contributions from the interstellar ra- 
diation field (ISRF) to the temperature distribution at the 
edge of the cloud. Figure 9 shows that To{r) is approx- 
imately a power law. At small radii, T£)(r) deviates by 
rising more steeply towards the center than a strict power 
law. Additionally, the ISRF causes an upturn in To (r) at 
large radii (about 1 pc for a 10^ L© source). The ISRF 
makes a siibstantial contribution to the temperature pro- 
file in regions forming low mass stars where the internal 
luminosity is low (Shirley et al. 2002a, Young et al. 2002) 
but has little effect in regions forming massive stars be- 
cause the temperature profile is dominated by the embed- 
ded source. While these sources may exist in regions of 
enhanced ISRF, the effects on the models are negligible 
even in our least luminous sources (L&oZ — 10^ Lq) unless 
the ISRF is a factor of 10 stronger than the standard value 
(see Evans et al. 2001 for a plot of the ISRF). For a more 
typical source with Lioi = 10^ Lq, the ISRF field must be 
at least 200 times stronger to change the best fit p by 0.25. 

For each of our models, the input density, n/, was nor- 
malized so that the model flux at 350 fim matched our 
observations, given a dust opacity at 350 fim. Dust opaci- 
ties were adopted from column 5 of the table in Ossenkopf 
and Henning (1994; hereafter 0H5), which were calculated 
for coagulated dust grains with ice mantles. These dust 



opacities (0II5) have been previously shown to match ob- 
servations of massive star formation regions by van der Tak 
et al. (1999, 2000) as well as low mass star forming regions 
(Evans et al. 2001, Shirley et al. 2002a, Young et al. 2002). 
We also considered opacities for grains without ice mantles 
taken from column 2 of Ossenkopf & Henning (1994; 0H2) 
because temperatures in regions with young stars may be 
high enough to destroy the mantles. The 0H2 opacities 
are higher at A > 350 fim {koh2/ nom = 1.8 at 700 /xm), 
but lower at shorter wavelengths. The crossover point is 
near 350 /xm, with koh2/koh5 = 1-1 at 350 /zm (see Fig. 
2 of Evans et al. 2001). Using 0II2 opacities thus results in 
best fit Uf lower by about 10%. Using 0II2 opacities also 
resulted in higher model fiuxes at long wavelengths, but 
lower fluxes at short wavelengths. No difference was found 
in the best fit p between the 0H2 and 0H5 opacities. 

For each value of p and n/, other model input parame- 
ters were adjusted to fit the observed values. The temper- 
ature of the star was taken to be the value corresponding 
to a star with L = L^oi (Thompson 1984). However, the 
results are very insensitive to the stellar temperature (van 
der Tak et al. 2000). The outer, Tq, and inner, rj, radii 
for each source are based on the angular extent of the ob- 
served radial profile and the chopper throw. We set To to 
be the sum of the chopper throw and twice the extent of 
the observed profile. In particular, ro is large enough to 
allow simulation of chopping. The inner model radius was 
taken to be Vq divided by 1000, so that it is small enough 
to be unresolved. The model's sensitivity to the radii was 
tested in M8E (see §4.1 below). The best fit values were 
found to be insensitive to both inner and outer radii (see 
§4.1, Shirley et al. 2002a, Young et al. 2002). 

The observation simulation program uses the tempera- 
ture distribution output from the radiative transfer code 
for the density model being tested. The code calculates 
observed fluxes and luminosities and generates a radial 
profile of normalized intensity. The model is convolved 
with the observed beam and chopping is simulated to pro- 
duce a more realistic radial profile for comparison with 
observations. For our models, an average one-dimensional 
representation of the actual beam was used (Fig. 2). The 
fit of the models with the observations was quantified by 
calculation of the reduced chi squared (Xr)- The Xr value 
for the radial profile follows the definition in Evans et al. 
(2001) and is denoted x^^a- We also compute a Xr value 
for the fit to the SED, denoted xIsed- ^he shortest wave- 
length point was generally left out of the calculation of 
XsED' because the model underestimates the flux at A < 
30 fim. This effect is well known in spherical models that 
do not account for holes and inhomogeneities in the cloud 
that allow these wavelengths to escape (e.g., van der Tak 
et al. 2000). 

The modeling scheme and the dependence of the de- 
rived quantities on the model parameters are discussed by 
Evans et al. (2002) and the sensitivity of the best fit p 
to uncertainties in other parameters is quantified for low- 
luminosity regions by Shirley et al. (2002a). They found 
that the largest source of uncertainty in p is the strength 
of the interstellar radiation field. For the hmiinous sources 
studied here, this is a minor source of uncertainty, as al- 
ready noted. 
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4.1. M8E: A Model in Detail 

M8E was used as a test case for checking the effects of 
the input parameters on the models. The 350 /im contour 
map of M8E is shown in Figure 4. A spectro-photometric 
distance of 1.8 kpc (Blitz, Fich & Stark 1982) was used. 
The observed bolomctric luminosity, Lboi, was 1.47 x 10* 
Lq. The best fit model yielded L^oi = 1-45 x 10* Lq, 
calculated from convolving the model emission with the 
beams used for the observations. The sensitivity of the 
model to the input internal luminosity was tested by de- 
creasing and increasing the parameter by a factor of two. 
While changing the luminosity affected the model SED by 
increasing or decreasing the flux at certain wavelengths, it 
had no effect on the radial profile fit. 

The fiducial density, n/, was fixed by matching the 
model flux at 350 /im to the observed value for each mod- 
eled density distribution. The best fit n/ = 1.2 x 10 
cm-3 for M8E. For M8E, = 4.3 x 10^ AU (2.1 pc) 
and rj = 430 AU. To test the sensitivity of the models 
to these parameters, ri was doubled and, independently. 
To was halved. Both tests yielded the same best fit p as 
the model with and Tq calculated in the standard way. 
Other resulting model values, such as L^oi and T^oi, also 
did not vary significantly. Therefore, we concluded that 
the models are insensitive to changes of a factor of two in 
r, and Tq. 

Figure 10 shows the best fit model for this source, a den- 
sity power law with p = 1.75 (xsed = 13.5, X350 = 0.52). 
As with a few of our sample, the best fit density distri- 
bution to the radial profile did not agree with the best 
fit to the SED for M8E. The Xsed ^as lower (8.3) with 
p = 1.5 and a lower fiducial density and central temper- 
ature. However, the slope of the model radial profile is 
clearly too shallow (x^c^a = 4.3, shown as a dashed line in 
Fig. 10). With M8E and the other sources that showed 
a disagreement in the best fit of the SED and radial pro- 
file, we report the best fit p to the radial profile because 
the shape of the opacity law and density inhomogeneities 
strongly affect the SED at wavelengths where the dust is 
optically thick. The radial profiles are much more effective 
in constraining the value of p. The p = 2.0 model (dotted 
line in Figure 10) produced a slope very similar to what 
was observed, but was not a best fit because the radial 
profile steepens at a smaller radius than observed. 

These models illustrate (Fig. 10) a conservative uncer- 
tainty in p of ±0.25; while p = 1.5 or 2 clearly do not fit 
the data as well as p = 1.75, they do lie at the edges of 
the error bars in the data. The X350 for models with p 
= 1.6-1.9 were near or less than one. The models were 
quantitatively distinguishable for 5p = ±0.2 {X350 = 2.4 
for p = 1.55). The limiting uncertainty in p appears to 
be the signal-to-noise of the profiles, rather than any sys- 
tematic effect. However, this uncertainty assumes that 
a power law density profile is a reasonable model; many 
sources are known to be more complex on smaller scales. 

Another source of uncertainty is the knowledge of the 
beam profile. Shirley et al. (2002a) tested the effects of 
the James Clerk Maxwell Telescope (JCMT) beam uncer- 
tainties and found 5p = Q.l. As were all sources, M8E 
was modeled with the average of the two observed beam 
profiles shown in Figure 2. For M8E, we also tried each 
of the two beam profiles, which were observed on different 



nights. Using the individual beam profiles did not change 
the result {6p < 0.05) for M8E. For these data, the beam 
uncertainties are negligible compared with those in the 
data, as long as the measured beam is used. If the beam 
is represented by a Gaussian, the effects are much larger 
(§6.1). 

4.2. Model Results 

Of our sample of 51 regions forming massive stars, 31 
could be modeled. We required that the source was not 
confused with multiple sources, that the map had high 
signal-to-noise, and that a range of flux density measure- 
ments were available in the literature. Most of our mod- 
eled sources arc well resolved, i.e., the deconvolved source 
size is at least the beam size, 9dec/(^mb > 1- We also 
consider some smaller sources to be resolved and require 
(^dec/Smb > 0.8 (or a source with an observed FWHM of 
1.3 times the FWHM of the beam) for modeling. Table 3 
lists 9dec/0mb for the sample. Three modeled sources have 
0.8 < OdeclOmb < 1 (W28A2, W43S, G31.41±0.31). The 
model results for these sources show steep density profiles 
(p = 2.25-2.5) and should be considered less robust than 
the rest of the sample, because they may not be resolved 
due to uncertainties in the beam. However, the steeper 
profiles may also be real; in §6.3 we discuss the correla- 
tion between p and Odec/^mb- Because of this correlation, 
our models arc biased against sources with p > 2.5 which 
would not be considered resolved. The requirement on 6dec 
also produces a bias against small sources. For a source at 
the median distance of 2.8 kpc, sources with rdec < 0.14 
pc would not qualify for modeling. If the minimum 350 
flux density of our modeled sample and the median 
distance are assumed, the minimum required mass for our 
models is 61 Mq. The sources with small angular extents 
were also in some cases the most luminous and distant 
sources (e.g., G12.21-0.10 with Lboi = 5.5 x 10^ Lg and 
D = 13.7 kpc). 

Figures 11 and 12 show the best flt model with the ob- 
served SED and radial profile for a sub-sample of sources 
that illustrate the full range of best fitting p values. Figure 
13 shows a histogram of the best fit p values. The mean 
and standard deviation of the histogram are (p) = 1.8±0.4. 
The standard deviation is about twice our estimated un- 
certainties on a single fit, providing marginal evidence for 
a range of actual values of p. Figure 13 also shows the 
distribution of densities at 1000 AU, n/. For this sample, 
(n/) = 1.2 x 10^ cm~^ and the median value is 1.4 x 10'' 
cm~^. In Section 6.3, we compare the distributions of p 
and n/ to those found in studies of low mass star forma- 
tion. 

The models generally flt the observed radial proflles 
very well. The value of X350 for the best flt models was 
less than one for nearly all the models. The majority of 
the x'sED values were under 10 over the range of wave- 
lengths where they were computed. For reasons described 
in §4, the model typically underestimates the emission at 
shorter wavelengths, where the dust is becoming optically 
thick. Models with density inhomogeneities (clumps, cavi- 
ties, etc.) or flattened structure might match the emission 
at shorter wavelengths, but such models introduce many 
free parameters. 

Recent studies of low mass star forming cores using the 
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same modeling techniques as presented here have found 
a correlation between p and the aspect ratio of the core 
(Shirley et al. 2002a and Young et al. 2002). Because our 
source shapes were affected by chopping, we cannot mea- 
sure reliable aspect ratios in enough sources to make this 
comparison. 

The modeled density power law exponent can also be af- 
fected by the inclusion of a compact central source, such as 
an UCH II region. For low mass star forming cores, the in- 
clusion of a compact source, in that case a disk, decreased 
p for the envelope by 0.5 (Young et al. 2002). 

The UCH II region at the position of the water maser 
in W3(0H) has a 3 mm flux density of 3.5 Jy (Wilncr ct 
al. 1995). For a conservative upper limit on the amount 
of flux the UCH II region contributes at 350 /xm, we have 
assumed that all of the 3 mm flux is due to ionized gas 
rather than dust and that the UCH II region is optically 
thick at 350 /im so that 5*,^ oc i^^. These assumptions give 
the maximum contribution of an UCH II region at 350 /im 
as 23% of the total observed flux. Including a compact 
source with this maximum contribution in our dust model, 
assuming that all of the emission from the UCH II region 
is included only in the central beam of our observations 
(as in Young et al. 2002), steepens the radial profile and, 
therefore, steepens the best fit p to the radial profile. For 
W3(0H), p increased by 0.3. The result is that the possible 
presence of a compact source introduces uncertainty into 
the models, because the modeled density power law of the 
envelope is steeper than the density structure if there were 
no compact source. We conclude that a UCH II region in- 
troduces uncertainty in our model results for the density 
structure of the envelope, 5p = —0.3, if it contributes more 
than 20% of the total 350 /im flux. 

Rick Forster at the Berkely-Illinois-Maryland Array 
(BIMA) generously provided recent 3 mm flux densities 
for several UCH II regions near the center of cores in this 
study. Three of the four sources have negligible contribu- 
tion from the UCH II region at 350 /zm. Therefore, the 
simple presence of a UCH II region does not necessarily 
imply an uncertainty in p. 

Our results for the density structure of massive star 
forming regions can be compared with theoretical predic- 
tions. The mean value of p (1.8) is incompatible with a 
logatropic sphere (p = 1; McLaughlin & Pudritz 1997) at 
about a 2 (7 level. However, both p = 1.5 and p = 2 are 
possible, so the inside-out collapse model of star formation 
(Shu 1977) cannot be ruled out. While a changing opacity 
as a function of radius might aff'ect the values of p and this 
conclusion, the opacity of a grain with (0H5) and with- 
out (0H2) an ice mantle is nearly the same at 350 lira, so 
evaporation of mantles will not have a large effect (see §4). 

5. MASSES 

5.1. Integrated Mass 

In a power-law mass distribution, the mass can only 

be defined within some specified radius. For sources with 
models, we can compute the integrated mass within a par- 
ticular radius, R, from 



Mint = 4:niJ.mnfr^- J r'^ ^dr 



= 47r/xmn/rP— ^j— p<3 , (4) 

where fi = 2.3, to is the hydrogen atom mass, is the 
gas density at r/ (1000 AU), and p is the best fit density 
power law exponent. These masses refer to the total mass 
of gas and dust within R. The mass is proportional to 
R^~P; for the mean value of p, Mint oc R^'^. 

One choice for R is rdec, the radius corresponding to the 
deconvolved source size. While this size has no intrinsic 
meaning if the density distribution is truly a power law, 
it is a fiducial size that is model independent. For the 
modeled sample, the integrated masses (Table 4) within 
rdec range from 6 to 1500 Mq and (r'Xjf) = 0.16 ± 0.09 
pc, the same as for the complete sample; (M(< rdec)) = 
250 ± 380 Mq with a median value of 120 M©, and 
{log(Af(< rdec))) = 2.0 ± 0.6 for this subsample. 

A second, more physical choice for R is the radius in- 
side which the density is actually enhanced over the sur- 
rounding cloud. This choice would define the total mass 
of the actual core, but we generally lack information for 
particular sources on the ambient density. In a study of 
extended cloud conditions in regions of massive star for- 
mation, including some studied here, AUers et al. (2002) 
find a typical ambient gas density of n 10^ cm^"^. The 
models for cores allow determination of r„ , defined to be 
the radius at which n = lO"' cm^'^. For the modeled sam- 
ple of 31 sources, (r„) = 0.42 pc, about 2.5 times {rdec)- 
Setting R = rn yields a core mass (M(< r„)) given in Ta- 
ble 4 for each core. Averaged over the cores with models, 
(M(< r„)) = 720 ± 860 M© and (log(M(< r„))) = 2.5 ± 
0.6. Figure 14 shows the distribution of sizes {rdec and r„) 
and integrated masses {M{< rdec) and M(< r„)) for the 
modeled sample. 

5.2. Isothermal Mass 

In order to estimate the mass for those sources without 
models, and hence a best fit p, the measured flux density 
was used to calculate a mass by assuming a single repre- 
sentative temperature. The isothermal total mass, M^so, 
was calculated according to the equation: 



= 5.09 X lO-^M© 5,(Jy)£>2(pc)(e^i-^/^*- - 1)(5) 

which assumes a single dust temperature, Tiso- is the 
observed flux at 350 y^m in a 120" aperture and D is the 
distance. We assumed the 0H5 opacity, = 0.10 cm^ g~^ 
of gas and dust at 350 /xm. Using the modeled sources with 
integrated masses. Mint, within a 120" aperture, we can 
calculate what assumption about dust temperature would 
give the best agreement with the masses from the mod- 
els. The "isothermal temperature", Tiso, is given by the 
equation: 

rp ^ hu/k 

ln(l + {2hv^K,Mint{MQ))/{c^S,{iy)D-'{^c))) 

(6) 

(Shirley et al. 2002a). For all the other cores, we assumed 
the mean temperature derived from the modeled sources, 
(Tiso) = 29 ± 9 K. This method allowed us to estimate the 
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masses for the complete sample (see Table 4 and Figure 8). 
We find {Mi,o) = 2020 ± 4410 Mq, a median M,so = 400 
M©, and (log(Miso)) = 2.8±0.7. The distribution of Mj^o 
is skewed to lower masses, but it has a tail of very massive 
cores. 

5.3. Mean Column Density 

The mean column density, S, has been used by McKee & 
Tan (2002) (and Tan & McKec 2002) to derive fundamen- 
tal physical quantities in regions forming massive stars. 
The observed S is a key parameter in the determination 
of the mean pressure, accretion rate, and star formation 
time (McKee & Tan 2002). Tan & McKee (2002) use the 
virial masses and radii from Plume et al. (1997) and report 
Tiyir ~ 1 g cm~^. For comparison, we calculate E using 
the masses derived from the dust continuum for our sub- 
sample of Plume et al. (1997). For the modeled sample, 
using M(< r„) and r„, ^mod = 0.19±0.12 g cm~^ which is 
significantly lower than Tj^ir (Tan & McKce 2002). How- 
ever, expanding the calculation to the complete sample 
results in a more consistent value. Using the isothermal 
mass, Miso and r„ (or (r„) for the sources that were not 
modeled), E = 0.73 ± 1.7 g cm~^. The mass accretion 
rate, rri*. varies as Ti^/'^ and the star formation time, f*/, 
as Ti~^/'^; therefore, a lower E decreases the rn* and in- 
creases the t*/ given by McKee & Tan (2002). 

It is important to note the inverse dependence of the 
masses and the mean column density on the opacity. 
Shirley et al. (2002b) compare the masses determined from 
the density distributions reported here with virial masses 
and find that the virial masses are larger on average by 
3.4. While this is quite good agreement considering that 
the opacities from different dust models can vary up to a 
factor of 10 (Ossenkopf & Henning 1994), the masses and 
mean column density may be a factor of three larger. 

6. DISCUSSION 

6.1. Comparison with Previous Models 

The average density distribution exponent, 1.8, is sig- 
nificantly higher than the 1.0 to 1.5 reported by van der 
Tak et al. (2000). Steeper best fit power laws were also 
found for all of the sources in common. The van der Tak 
et al. (2000) best fit models were based on molecular line 
emission observations, and were not necessarily the best 
fits to their dust continuum emission (see their Figure 8). 
However, there is still a discrepancy between their mod- 
eled 350 /im radial profiles and our models even though 
the same method was used. The difference stems from 
convolving the model with an observed beam rather than 
a 10" Gaussian, as was done with the van der Tak et al. 
(2000) models. Figure 15 compares a modeled radial pro- 
file for GL2591 using a 10" Gaussian and the observed 
beam with p = 2.0 and all other parameters the same. The 
Gaussian beam requires a much shallower density power 
law to fit the observations. In the case of GL2591, van der 
Tak et al. (2000) reported a best fit based on molecular 
line emission oi p = 1.0 and show that the dust emission 
radial profile would be fit by p ~ 1.25 (see their Figure 8). 
We find p = 2.0; this was the largest discrepancy caused by 
the beam for the overlapping sample. On average, van der 
Tak et al. (2000) report a shallower power law by 0.4. If a 



14" Gaussian, which we find better describes the beam at 
the time of our observations, is used, the difference is also 
5p = —0.75 for GL2591; the main effect is in the sidelobes. 

6.2. Comparison with Other Studies 

6.2.1. Hatchell et al. 2000 

Hatchell et al. (2000) modeled the 450 and 850 /^m 
dust emission for five massive star forming regions se- 
lected by emission from dense molecular gas. Three of 
their sources are included in this study (G12.21— 0.10, 
G31.41-f0.31 and G13.87+0.28). Hatchell et al. (2000) 
fit the SED and radial profile at 450 and 850 /um using 
the dust radiative transfer code DUSTY (Ivezic & Elitzur 
1997). For G13. 87-1-0.28, the only modeled source in com- 
mon, Hatchell et al. (2000) report p = 1.5, which is shal- 
lower than our best fit p = 1.75. However, although p 
= 1.0 and 2.0 were ruled out by Hatchell et al. (2000), 
intermediate p values were not investigated in their study. 

Hatchell et al. (2000) find substantially higher masses 
(by a few times ten) than in our study. This discrepancy 
is due to the use of different dust properties. Hatchell et 
al. (2000) used Draine & Lee (1984) ice-free graphite and 
silicon dust grains as opposed to 0H5 coagulated grains 
with ice mantles. The Draine & Lee (1984) opacities are 
much smaller and result in higher mass estimates. 

Hatchell et al. (2000) included compact central cores in 
some of their models to fit the radial profiles. They found 
that p = 1.5 density distribution was too shallow, but a 
central core with & p = 1.5 envelope fit the data. Hatchell 
et al. (2000) used the core model for G12.21-0.10 and 
G31.41-I-0.31, because the 450 and 850 /um radial profiles 
fell off steeply within 20". The 350 /im radial intensity 
profiles for these sources are also steeper than average. 
We report p = 2.25 for G31. 414-0. 31 although it was not 
well resolved with Odec/Omb = 0.8. G12.21-0.10 did not 
fit our modeling criteria because 6dec/&mb < 0.5, so we 
considered it unresolved. 

6.2.2. Beuther et al. 2002 

Beuther et al. (2002) studied the density structure of 69 
regions forming massive stars (Sridharan at al. 2002) with 
1.2 mm continuum and CS emission. Their resolution at 
1.2 mm was very similar to ours at 350 /xm, providing an 
interesting comparison. To determine the density struc- 
ture, they fit the radial intensity profiles with a broken 
power law, steeper in the outer region, and assumed a 
power law temperature distribution (T oc r^*^'^) to deter- 
mine the density power law. We tested the effect of the 
power law temperature distribution by modeling a source 
with T oc r^" ** rather than the temperature distribution 
calculated from the dust code (see §4). Assuming a power 
law resulted in less flux at shorter wavelengths, up to 60% 
at 25 /im, and a 20% higher n/. However, the radial profile 
was unaffected and the best fit p did not change. There- 
fore, the resultant density distributions from these studies 
can be reasonably compared. 

Beuther et al. (2002) fit the radial intensity profiles on 
average with / oc r~^ '^ within 32" and r~^-^ in the outer 
regions. They report a mean density power law index de- 
rived from the mean inner radial index of {p) = 1.6 ± 0.5. 
These results are consistent with our findings, although 
we have modeled the core envelopes well beyond 32". As 
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noted in §3, we found no systematic tendency for intensity 
profiles to steepen beyond 40". This difference between the 
two studies could reflect differences in the sample or differ- 
ences in the observational details. For example, our chop 
throw was about twice that used by Bcuthcr ct al. (2002). 
More fundamentally, the agreement in the inner regions is 
very reassuring, indicating that observations from 350 to 
1200 /um are tracing the same distribution in these sources. 

Their sample has a mean mass, based on the integrated 
flux density and using either the unambiguous or near dis- 
tance, of (M) = 1550 Mq. This should be compared to 
our {Miso) = 2020 Mq, which was derived in a similar way. 
However, they used different assumptions about opacities. 
Beuther et al. (2002) explain that the use of 0H5 opacities, 
as used in our study, "would result in masses and column 
densities about a factor 4 lower". The factor of 4 makes 
their average masses about one fifth of our mean value or 
about equal to the median Miso (397 Mq). The Bcuthcr 
et al. (2002) "0H5 scale" mass distribution is shifted to 
lower masses than the distribution of isothermal masses in 
this study (sec Fig. 8) and peaks near log(M) = 2.2 com- 
pared to log(Miso) = 2.8. The discrepancy in the masses 
can further be explained by the use of different apertures 
in the measurement of flux densities between the two sam- 
ples. Beuther et al. (2002) measure flux densities within 
the 5% level of the peak emission. In our study, the aper- 
ture was generally larger, 120", resulting in larger flux 
densities. Since the isothermal mass is proportional to the 
flux density, we expect higher masses than Beuther et al. 
(2002). 

Different methods were also used to calculate the lumi- 
nosities. Sridharan et al. (2002) calculate the luminosities 
of the same sample studied by Beuther et al. (2002) by in- 
tegrating a two-component greybody curve fit to Infrared 
Astronomical Satellite {IRAS) data. They find an aver- 
age luminosity {{L) = lO** Lq) a factor of 25 lower than 
the mean luminosity reported here from integrating the 
observed SEDs, which in many cases also included IRAS 
fluxes. However, only a negligible fraction of this discrep- 
ancy is the result of different methods. For the two sources 
in common, Sridharan et al. (2002) find about a 20% lower 
luminosity than Lobs when the sources are placed at the 
same distances. Because of their selection against UCH II 
regions, the Sridharan et al. (2002) sample may contain 
younger and, therefore, less luminous sources than our 
sample. Sridharan et al. (2002) report a lower luminosity 
to mass ratio in their sample than that found in UCH II 
regions, implying that the ratio increases as a core evolves. 
Section 6.4 discusses the luminosity to mass ratio of our 
sample. 

6.3. Comparison to Low Mass Star Formation, 
Evolution, and Correlations 

Figure 13 shows a histogram of p values for low mass star 
forming regions (Young et al. 2002 and Shirley et al. 2002a) 
alongside the distribution of p for the massive stars in this 
study. Both samples were modeled with the techniques 
described here. For the low mass cores, {p) = 1.6 ± 0.4 
and increases to (p) = 1.8 if cores with high aspect ratios 
are left out of the mean (Young et al. 2002). Although the 
modes of low and high mass star formation are drastically 
different, the distributions of p are strikingly similar. The 



mean fiducial density in the regions forming massive stars 
((n/) = 1.2 X 10* cm~'^) is over two orders of magnitude 
greater than that for the low mass cores {{nf)iowmass = 
5 X 10^ cm~^) (Figure 13). Because the fiducial density is 
referred to 1000 AU, which we do not resolve, the absolute 
value should not be taken too literally; however, the sim- 
ilar values of p imply that the densities are higher at all 
radii by about two orders of magnitude in these sources 
than in cores forming low mass stars. The similarity in p 
suggests that similar processes determine the shape of the 
density distribution in regions forming low and high mass 
stars. 

In order to address the problem of a coherent evolution- 
ary sequence for regions forming massive stars, we exam- 
ine quantities commonly used as indicators of evolution 
in regions forming low mass stars. One such indicator is 
the bolometric temperature, T^oi (Myers & Ladd 1993), 
the temperature of a blackbody with the same mean fre- 
quency as the observed SED. The range of bolometric tem- 
peratures in this study was relatively small, 46 - 173 K, 
compared to that found by Myers & Ladd (1993) for low 
mass objects, almost two orders of magnitude. Another 
indicator is the ratio of total to submillimetcr luminos- 
ity, Lf,oi/ Lgmm- For low mass star forming regions, both 
quantities are used as evolutionary indicators to mark the 
boundary between Class and Class I sources; Ltoi /Lsmm 
increases as a source evolves, and Andre et al. (1993) de- 
scribed sources with Lj^oi/Lsiyim < 200 as Class 0. Chen 
et al. (1995) defined Class objects to be sources with 
Tboi < 70 K. Tboi is plotted against LhoilLsmm in Figure 
16, which shows that while the Class and I definitions 
may not translate directly to massive stars, there is some 
correlation between the two measures. The trend implies 
that high bolometric temperature may be an indicator of 
more evolved sources. Young et al. (2002) plot Thoi versus 
I^boi I Lsrnrn for both low mass sources and this sample and 
find that the Class O/I boundary is not as clearly defined 
in regions forming low mass stars as in high mass regions. 

Since Thoi may be an indicator of evolution, we looked 
for correlations with other parameters. Figure 16 shows no 
significant correlation between Ti,oi and p or Uf. Another 
indicator of evolution is the ratio of envelope to stellar 
mass {Miso /Mi,). As a source evolves and more mate- 
rial is accreted onto the central source, the ratio should 
decrease. To approximate the stellar mass, M*, we as- 
sumed that the luminosity is dominated by a single mas- 
sive star, and that oc L^/^-^. This mass- luminosity 
relationship is for zero age main sequence stars and has 
been used for a range of stellar masses (van der Tak et al. 
2000, Shirley et al. 2002a). Observations of binary stel- 
lar masses and luminosities suggest that this relationship 
is valid to approximately 60 Mq (Scalo 1986). The lumi- 
nosities in this sample suggest very few of these cores have 
a central source larger than 60 Mq, and the sources with 
the highest luminosities are likely forming more than one 
massive star. Figure 16 also shows that for regions forming 
massive stars there is not a significant correlation between 
Tboi and Miso/M^,. However, the high values for T^oi (> 
80) occur for MisojMt, < 70, indicating a more evolved 
source. 

Figure 17 plots Miso/M^, and another potential evolu- 
tionary indicator, the far infrared color (Feo/Fioo), which 



Physical Conditions for Massive Star Formation 



9 



is the ratio of flux densities at 60 /im and 100 /xm, both 
observed with the same instrument [either IRAS or Kuiper 
Airborne Observatory (D. Jaffe 2001, private communica- 
tion)]. Our data do not show the same direct correlation 
between far-infrared color and ratio of envelope to stel- 
lar mass as reported by van der Tak et al. (2000). How- 
ever, their conclusion that bluer far-infrared colors only 
occur for lower mass ratios does apply. A more evolved 
source, as measured by an envelope to stellar mass ra- 
tio less than 70, may have a T^oi > 100 and far-infrared 
color, Fqq/Fxqo > 0.7, while these values are never found 
for Miso/M^ > 70. Figure 17 also plots p versus Miso/M^ 
and indicates that less evolved sources (i.e., high Miso/M^,) 
may have steeper than average density distributions. 

Possible correlations with the density distribution were 
also examined. Van der Tak et al. (2000) found no correla- 
tions between envelope mass or the internal luminosity and 
the density law exponent, p. However, Figure 18 shows 
our larger sample does suggest some correlation with each 
of these parameters (r ~ 0.5). The density distribution 
steepens with increasing mass and luminosity. Figure 18 
also plots the ratio of the deconvolved source diameter to 
the FWHM beam size, 9dec/()mb, and the distance versus 
p. There is a strong relationship between Odec/Omb and the 
best fit p. Better resolved sources tend to have shallower 
density profiles, as do more nearby sources. The correla- 
tion of Odec/Omb with p could bc the result of a systematic 
effect in the models or resolution. Three of the six sources 
with p> 2 (GlO.60-0.40, G31.41+0.31, and W43S) are 
also the most distant, suggesting that better resolution is 
needed to confirm their density distributions. The cor- 
relation of 9dec/dmb with p could also indicate a physical 
difference in the cores, such as different initial conditions. 
Young et al. (2002) found the same correlation for low 
mass cores where distance is not the same limiting factor 
as in this study. A source with a steeper density distribu- 
tion (higher p) naturally has a smaller size {Odec) and will 
be less well resolved (Figure 18(c)). 

6.4. Luminosity to Mass Ratios 

6.4.1. Comparison to Galactic Studies 

The luminosity to mass ratio, L/M , is often used as a 
measure of the star formation rate per unit mass. Most 
studies in our Galaxy have used the mass of the entire 
molecular cloud, determined from CO emission. The re- 
sulting average [L/M = 0.4 Lq/M©) is quite low (Bronf- 
man et al. 2000). For molecular clouds that contain H II 
regions, the luminosity to mass ratio from CO emission is 
an order of magnitude higher, L/M = 4 Lq/Mq (Mooney 
& Solomon 1988), but the dispersion is over 2-3 orders of 
magnitude (Evans 1991; Mead, Kutner, & Evans 1990). 
CO emission traces the less dense gas of the entire molec- 
ular cloud, while the dust continuum emission at submil- 
limeter wavelengths measures the mass of very dense gas 
actually involved in star formation. We computed val- 
ues of Lboi/Miso because Miso could be obtained for the 
largest number of sources. Figure 19 plots log{Lboi / Miso) 
versus log(Miso). Both L^oi and Miso are proportional to 
the square of the distance, so the distance uncertainties 
are minimized in the ratio. The ratio and, hence, the star 
formation rate per unit mass is constant over the entire 
mass range. The dispersion of log{Li,oi / Mi^o) is about one 



and a half orders of magnitude, which is significantly less 
than the 2-3 orders of magnitude when the mass is traced 
by CO. 

Figure 19 also plots \og{Lhoi) versus \og{Miso) and shows 
that the luminosity to mass ratio traced in massive cores 
is significantly higher than the average L/M in molecular 
clouds with H II regions. Figure 20 shows the distribution 
of \og{Lboi/Miso), which peaks near (log (Lbo/Z-^iso)) = 
2.0 ±0.4; the median log( Lboi / Miso) is also 2.0. The mean 
value of the ratio is {Ltoi/M^so) = 140±100 Lq/Mq, with 
a median value of 120 Lq/Mq. This {L^oi/Miso) is 30 
times that derived from CO for Galactic molecular clouds 
containing H II regions. Clearly the star formation rate 
per unit mass of dense gas is much higher, indicating that 
the dense gas traced by submilhmeter emission is the lo- 
cation of massive star formation. 

The histogram of L^oi / Miso in Figure 20 shows that the 
distribution is strongly skewed. Most sources have low ra- 
tios, but there is a tail of very high ratios, reaching up to 
490 L0/M0. A few sources could be affected by confusion 
with multiple sources in the large IRAS beam, resulting 
in a luminosity that is too high. This problem is difficult 
to avoid since IRAS is often the only source of mid- to 
far-infrared fluxes. Higher spatial resolution is needed to 
eliminate multiple source confusion in flux measurements. 
The Stratospheric Observatory For Infrared Astronomy 
(SOFIA) that will fly later in the decade will be ideally 
suited to address this problem and fix the upper limit to 
L/M for Galactic sources. 

The mean L/M for the Beuther et al. (2002) sam- 
ple {{L/M) = 20 ± 18 L0/M0, Sridharan et al. 2002) 
is 7 times lower than in this study {{L^oi/ Miso) — 140 
Lq/Mq). If we decrease their masses by a factor of four, 
putting them on the "0H5 scale", the values are in bet- 
ter agreement {{L/M) ~ 80 Lq/Mq). For sources with a 
distance ambiguity, we assume their near distance, which 
may produce a downward bias in luminosities and masses. 
Eliminating this bias by averaging only sources in their 
sample for which the distance is not ambiguous, gives a lu- 
minosity to mass ratio {{L/M) = 120 ±90 L0/M0) that is 
consistent with the results of our study. However, Sridha- 
ran et al. (2002) report an L/M significantly smaller than 
that of a sample of UCH II regions (Hunter 1997, Hunter 
et al. 2000). They suggest that the cores in their sample 
are in a younger pre-UCH II phase and that L /M increases 
as the cores evolve and develop UCH II regions (Sridharan 
et al. 2002). Hunter et al. (2000) describe their sample as 
having UCH II region far-infrared colors or thermal radio 
continuum emission. If we also convert the opacities used 
by Hunter et al. to the 0H5 scale, then the average lumi- 
nosity to mass ratio for their sample is {L/M) = 280 ± 370 
Lq/M© with a median of 150. The mean L/M for the 
UCH II regions studied by Hunter (1997) and Hunter et 
al. (2000) is higher, but the median is not inconsistent 
with the results of the work presented here for cores both 
with and without UCH II regions. 

In order to test the hypothesis that the luminosity to 
mass ratio increases as a massive star forming region de- 
velops an UCH II region, we divided our sample and calcu- 
lated L})oi / Miso for cores with and cores without UCH II 
regions. We do not find the same disparity in L/M be- 
tween our two subsamples as between the results of Srid- 
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haran ct al. and Hunter (1997 & Hunter et al. 2000). For 
the cores associated with UCH II regions, {L^oil Migo) = 
170 ± 130 (with a median of 120), which is lower than the 
mean of the sample of Hunter and Hunter et al., but within 
the dispersions of the two samples. The mean ratio of the 
remaining cores that do not contain an UCH II region is 
not significantly lower, {Lbd / Miso) = 130 ±90 with a me- 
dian of 100. However, in a CS study with a larger sample 
of which ours is a subset, Shirley et al. (2002b) found that 
the median luminosity to virial mass ratio is more than a 
factor of two higher in cores with UCH II regions than in 
cores without UCH II regions. 

6.4.2. Comparison with Extragalactic Studies 

Submillimeter wavelengths are becoming increasingly 

important in the the study of star formation in other galax- 
ies, especially at high redshifts (Blain et al. 2002). A recent 
submillimeter survey measiircd the dust masses, Mo. and 
far-infrared to submillimeter (~ 1 to 1000 /xm) luminosi- 
ties, LpiR, of bright local IRAS galaxies (Dunne & Bales 
2001). Dust is often used as a tracer of mass for more dis- 
tant sources as well (e.g., Omont et al. 2001, Calzetti et al. 
2000, Benford et al. 1999). These studies assume a single 
dust temperature, usually 50 K, based on the far-infrared 
and submillimeter SEDs of the galaxies (e.g, Calzetti et 
al. 2000, Benford ct al. 1999). The average T,so of 29 K 
found in this study is substantially lower. If we compared 
the mass within r^ec, we found that T^so was about 50 
K, but when using the full extent, much cooler dust con- 
tributed substantially. The dust mass, Md, of galaxies is 
usually calculated in the same way as our Miso (Equation 
(5)). However, different opacities are used and only the 
mass of the dust is calculated, not the mass of gas and 
dust as with Miso- 

For high redshift (z > 4) quasars, Omont et al. (2002) 
report a mean (and median) dust mass, {Mu) ~ 5 x 10^ 
Mq for a range of 2 x 10* to 1 x 10^ Mq. A typical 
far-infrared luminosity (> 50 /xm) for their sample is 
Lfir - 10" L©, giving Lfir/Md = 2 x lO" Lq/Mq. 
For purposes of comparison, we calculated an isothermal 
dust mass, M^, using the same temperature (50 K) and 
opacities as Omont et al. (2002). For the regions in this 
study {Lobs/Mo) = (1.4±1.0) x 10^ Lq/Mq ranging from 
2 X 10^ to 4.5 x 10'' L0/M0. The luminosity to mass ra- 
tio for high redshift quasars is similar to the higher values 
in our sample, suggesting that starbursts form stars as if 
most of their molecular material acts like the most extreme 
regions forming massive stars in the Milky Way. However, 
starbursts might not be the only source of L pj r in distant 
quasars. Omont et al. (2002) suggest central active galac- 
tic nuclei as another source of dust heating. In that case, 
the star-forming L/M would be smaller. 

7. SUMMARY 



We have presented dust continuum maps of 51 regions 
forming massive stars with a large range of sizes and 
masses, {raec) = 0.16 ± 0.10 pc and (Mj^o) = 2020 ± 4410 
Mq. We find that the peak of the dust emission is more 
often coincident with the water maser position than the 
UCH II region, implying that the dust may better trace 
the earlier stages of massive star formation. Modeling 
a subset of 31 sources yielded a mean density index for 
power laws (n(r) = (r/r/)-P) of (p) = 1.8 ± 0.4. The 
dispersion is about twice the expected uncertainty for an 
individual source, suggesting some real dispersion among 
sources. The mean value for p is incompatible with a loga- 
tropic sphere (p = 1), but p = 1.5 or 2 are possible. The 
mean value and variation in p are similar to those found 
for low mass regions (Young et al. 2002). 

For the modeled sources, integrated masses within two 
different fiducial radii were presented in addition to the 
isothermal masses. The mean mass within the 350 fim 
half-power radius (rdec) is (M(< rdec)) = 250 ± 380 Mq. 
For the larger, more physical radius where the density falls 
off to the ambient level (n = 10^ cm-^), (M(< r„)) = 720 
± 860 Mq. The mean column density for the complete 
sample was found to be E = 0.73 ± 1.7. However, because 
of the inverse dependency of mass on opacity, compari- 
son with virial masses suggests that the masses and mean 
surface density could be about a factor of three larger. 

The density structure of massive star-forming cores was 
found to be consistent with Beuther et al. (2002). The 
results were also consistent with van der Tak et al. (2000) 
when the difference in a Gaussian and the observed beam 
profile was taken into account. Our data also confirm some 
of the conclusions made by van der Tak et al. (2000) with 
respect to possible evolutionary indicators. There is a 
trend of rising T^oi with increasing L/Lgmm, suggesting 
that either of these quantities could be tracing evolution. 

The luminosity to mass ratio, a tracer of star formation 
rate per unit mass, has a mean of {Lboi/Miao) — 140 ± 
100 L0/M0 for dense gas traced by dust emission, with a 
tail extending up to about 500 Lq/Mq. These values are 
much higher than those based on masses from CO emis- 
sion. The mean luminosity to mass ratio derived using the 
dust mass is similar to that in extreme starburst galaxies, 
especially those seen at substantial redshift. 
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Fig. 1. — The distribution of distances in our sample. N is the number of sources. 
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Fig. 2. — Radial profiles of Uranus on 1998 July 23 and 1998 July 25 used to measure the beam plotted with a 14" Gaussian. The error 
bars represent the weighted mean variation in the annulus for which each point of the profile is calculated. We adopt a FWHM beam size, 

drub, of 14". 
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Fig. 3. — 350 dust emission contour maps. The arrows indicate tlie cliop direction but are not indicative of the chop length (~ 100"). 
The plus signs indicate the positions of known UCH II regions (W3(OH): Wilncr ct al. 1995). The contour levels are as follows: IRAS 
00338+6312 {4a), G123.07-6.31 (3(7, then in increments of 20% (lOtr) of the peak), W3 (4<t), W3(0H) {3a, 20% {Wa)), IRAS 02395+6244 
{2a), IRAS 02461+6147 {2a), G137.07-3.00 {2a), GL490 {3a, 20% {6a)), Ori-IRC2 {3a, 10% {7a)), S231 (4ct), S235 {2a), S241 {2a). 
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Fig. 4. — 350 fim dust emission contour maps. The arrows indicate tlie cliop direction but are not indicative of the chop length (~ 100"). 
The plus signs indicate the positions of known UCH II regions (MonR2, W28A2, G8.67, G10.60, G12.21: Wood & Churchwell 1989; RCW142: 
Walsh et al. 1998; G12.42: Jaflfc ct al. 1984). The contour levels arc as follows: MonR2 (4o-), S252A (4o-), RCW142 (Scr, then in increments 
of 20% (260-) of the peak), W28A2 (3(7, 20% (29(7)), M8E (3(7, 10% (5a)), G9. 62+0. 10 (3a, 20% (28(7)), G8.67-0.36 (3(7, 20% (25a)), 
GlO.60-0.40 (3(7, 20% (30a)), G12.42+0.50 {3a, 20% (17(7)), G12.89+0.49 (3a, 20% (12a)), G12.21-0.10 (3(7, 10% (20(7)), G13.87+0.28 (3a, 
10% (9a)). 
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Fig. 5. — 350 /xm dust emission contour maps. The arrows indicate the chop direction but are not indicative of the chop length (~ 100"). 
The plus signs indicate the positions of known UCH II regions (Wood & Churchwell 1989). The contour levels axe as follows: W33A (4it), 
G14.33-0.64 (3(7, then in increments of 20% (17ct) of the peak), GL2136 (3ct, 20% (8(t)), G19.61-0.23 {3(J, 20% (33o-)), G23.95+0.16 (3(t, 
20% {5a)), G24.49-0.04 {3a, 20% {15a)), W43S {3a, 20% (29ct)), G31.41+0.31 {3a, 20% {18a)), G40.50+2.54 {3a, 20% (IOct)), G35.58-0.03 
(4ct), G45.07+0.13 {3a, 20% {6a)), G48.61+0.02 {3a, 10% {5a)). 
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Fig. 6. — 350 jum dust emission contour maps. The arrows indicate the chop direction but are not indicative of the chop length (~ 100"). 
The plus signs indicate the positions of known UCH II regions (W51M: Scott 1978; S88B: Wood & Churchwell 1989; S87, ONI, S106, S157: 
Kurtz et al. 1994; GL2591, S140: Tofani et al. 1995). The contour levels are as follows: W51M (3<7, then in increments of 20% {19<t) of the 
peak), S87 (4o-), S88B (4(t), ONI {3a, 20% (Ua)), ON2S {4a), S106 {4a), GL2591 {3a, 20% (Hct)), G97.53+3.19 {2a), BFSllB {2a), S140 
(4ct), CEP A {3a, 10% {7a)), S157 {2a). 
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Fig. 7. — Composite 350 contour map of NGC 7538. The arrows indieate the ehop direction but are not indicative of the chop length 
(~ 100"). Clockwise from left are the sources IRS9, IRSl, and IRSll. The visible H II region NCG7538 (S158) hes north of IRSl. The plus 
sign indicates the position of the UCH II region associated with IRSl (23'' 11"" 36.6" 61° 11' 36.6", Wood & Churchwell 1989). The contour 
levels are as follows: IRS9 (4ct), IRSl (3(t, 10% (6(t)), IRSll (3(t, 10% {6a)). 
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Fig. 8. — Distribution of radii and masses for the complete sample, r^ec is the deconvolved half-power size as determined from the 350 /um 
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3 4 5 

Log r (AU) 

Fig. 9. — The solid line is the model dust temperature distribution for internal source with Lf,gi = 10** L©. The dashed line is a Tu{t) oc r~^ '^ 
power law. The modeled distribution deviates from a strict power law at large r due to the ISRF and at small r because of optical depth 
effects. 



Physical Conditions for Massive Star Formation 



21 




Log A (Atm) Log b (AU) 

Fig. 10. — The left figure shows the observed (circles with error bars) and best fit model (soUd line) SED (xfj = 14) for M8E. On the right, 
the solid line is the best fit model for the observed normalized radial intensity profile (error bars), p = 1.75 (solid line, xfj = 0.52). The error 
bars in the radial profile represent the weighted mean variation in the annulus for which each point of the profile is calculated. Also plotted 
are the p = 1.5 (dashed line, Xjj = 4.3) and p = 2.0 (dotted line, Xjf = 3-4) model profiles. 
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Fig. 11. — Observed SEDs (circles) and normalized radial profiles (error bars) for a subsamplc of sources witfi best fit models (solid line) 
for p = 0.75 to 1.75. The error bars in the radial profiles represent the weighted mean variation in the annulus for which each point of the 
profile is calculated. X%ed calculated for SED points > 12 //m. 
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Fig. 12. — Observed SEDs (circles) and normalized radial profiles (error bars) for a subsamplc of sources with best fit models (solid line) for 
p = 2.0 to 2.5. The error bars in the radial profiles represent the weighted mean variation in the annulus for which each point of the profile 
is calculated. Xsbd calculated for SED points > 12 //m. 
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Fig. 13. — Histograms of fiducial densities (n/) and best fit density power law exponents (p) for high (left) and low (right) mass star 
forming cores modeled using the same methods. For the massive cores, (rij-) = 1.2 X 10* em^'* and (p) = 1.8 ± 0.4. For the low mass cores, 
(n/) = 5.2 X 10^ cm-=^ and (p) = 1.6 ± 0.4; (p) = 1.8 if sources with high aspect ratios are not included (Young et al. 2002). 
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Fig. 14. — Distributions of radii and integrated masses for the modeled sample. The top two panels are the histograms of ''2^°'' and the 
mass within rdec- The bottom panels show the distributions of rn, the radius at which the density falls to 10* cm~^, and the mass within 
that radius. 
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Fig. 15. — Model radial profile with p = 2.0 using the observed beam profile (solid line) and a 10" Gaussian beam (dotted line) as used 
by van dcr Tak ct al. (2000) plotted with the observed radial profile of GL2591 (error bars). The error bars in the radial profiles represent 
the weighted mean variation in the ammlus for which each point of the profile is calculated. The use of a 10" Gaussian beam decreases the 
modeled p for GL2591 by 0.75, accounting for the discrepancy in the value of p reported by van der Tak et al. (2000) (fs 1.25) and in this 
study (2.0). 
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Fig. 16. — Relationships with Tijoi- (a) Relationship between bolometrie temperature and the ratio of total to submillimeter luminosity, 
both eonsidered indicators of evolution in low mass star forming cores. The solid lines indicate the divisions between Class and I low mass 
protostars for each indicator. Class I sources have T50; > 70 K and L/Lsrnm > 200. The linear correlation coefficient, r, is 0.52. (b) Tjoj 
versus another evolutionary indicator, Miso/M^ where M» = L^Z-^-^ (j- = 0.06). More evolved sources (Migo/Mt < 70) have higher Tj,o( (> 
80). (c) Tfigi versus the density power law index p (r = —0.25). (d) Tfi^i versus the fiducial density, ny, for the modeled sample (r = —0.33). 
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Fig. 17. — Relationships with the envelope to stellar mass ratio. Left: the far-infrared color versus Mjgo/M* (r = 0.07). More evolved 
sources (Mjgo/M* < 70) have bluer colors (> 0.7). Right: p versus Miso/M* {r = 0.31). 
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Fig. 18. — Relationships with p: (a) the density power law index p versus the logarithm of the mass within the 350 /im half-power radius 
(r = 0.51). (b) p versus the logarithm of the modeled internal luminosity (r = 0.49). (c) p versus the ratio of FWHM source size to beam 
size (r = —0.86). Only sources with ddec/Smb > 0.8 were modeled. Smaller sources have steeper density profiles. The solid line represents 
models of dust emission with p = 0.5 to 2.5. Sources significantly above the line (e.g., S140 and S88B) have extended aspherical emission, 
(d) p versus distance (r = 0.60). 
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Log M,,„ (MJ Log M,,, (M^) 

Fig. 19. — Left: Plot of log{Li,oi / Miso) versus log(Afiso). The dispersion in \og(Li,ci/ Migo) is about 1.5 orders of magnitude, significantly 
less than in CO observations (2-3 magnitudes). The error bar shows the full range from the CO observations (Evans 1991). Right: Plot of the 
logarithm of Lijoi versus the logarithm of M^so {i" = 0.89). The solid line is the least squares fit to the data log(Li,o;) = 1.9 + \og(Mii,o)- The 
dashed line is based on CO masses for molecular clouds containing H II regions, log(Lj,o() = 0.58 + log(Mjso) (Mooney & Solomon 1988). 
The luminosity to mass ratio is much higher in massive star forming cores than in molecular clouds. 
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Fig. 20. — Distribution of the luminosity to mass ratio, L},ai/ Miso, (left) and log {Li„ij Miso) (right). The histogram of Li„i/Mi, 
a tail of high ratios out to 490 /Mq . 
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Table 1 



Observed Properties 



Source 


a (1950.0) 


6 (1950.0) 

(° ' ") 


D 


Centroid'^ 


350/um 




D 






(kpc) 


("/') 


(Jy) 


{") 


Ref. 


IRAS 00338+6312 


00 33 53.3 


63 12 31 


0.85 


(-10,-5) 


160+32 


30 


10 










410+82 


120 




G123.07-6.31 


00 49 29.2 


56 17 36 


2.2 


(-15,-9) 


160+32 


30 


6 










290+58 


120 




W3^''' 


02 21 53.1 


61 52 20 


2.3 


(-30,9) 


30+6 


30 


7 










220+44 


120 




W3(0H) 


02 23 17.3 


61 38 58 


2.4 


(7,-2) 


400+80 
1130+230 


30 
120 


6 


IRAS 02395+6244 


02 39 31.0 


62 44 16 


8.1 


(-9,12) 


20+4 


30 


5 


IRAS 02461+6147 


02 46 11.7 


61 47 34 


4.5 


(-9,4) 


20+4 


30 


5 


G137.07-3.00 


02 54 11.2 


56 17 36 


4.9 


(1,10) 


20+4 


30 


6 


GL490'' 


03 23 38.9 


58 36 33 


0.9 


(-3,2) 


80+16 


30 


11 










180+36 


120 




Ori-IRC2 


05 32 47.0 


-05 24 24 


0.45 


(-4,3) 


1940+390 


30 


5 










7680+1540 


120 




S231 


05 35 51.3 


35 44 16 


2.3 


(0,0) 


190+38 


30 


6 










522+100 


120 




8235'= 


05 37 31.8 


35 40 18 


1.6 


(1,11) 


50+10 


30 


6 










240+48 


120 




S241 


06 00 40.9 


30 14 54 


4.7 


(8,9) 


30+6 


30 


6 










40+8 


120 




Mon R2'' 


06 05 17.0 


-06 22 40 


0.9 


(15,10) 


150+30 


30 


6 










1400+280 


120 




S252A 


06 05 36.5 


20 39 34 


1.5 


(2,6) 


130+26 


30 


6 










320+64 


120 




RCW142 


17 47 04.0 


-28 53 42 


2.0 


(-5,8) 


530+110 


30 


6 










670+130 


120 




W28A2^ 


17 57 26.8 


-24 03 54 


2.6 


(0,-3) 


950+190 


30 


6 










1580+320 


120 




M8E 


18 01 49.1 


-24 26 57 


1.8 


(-8,3) 


210+42 


30 


6 










380+76 


120 




G9.62+0.10 


18 03 16.0 


-20 32 01 


5.7 


(0,3) 


590+120 


30 


4 










1150+230 


120 




G8.67-0.36 


18 03 18.6 


-21 37 59 


4.5 


(1,4) 


650+130 


30 


8 










1160+230 


120 




GlO.60-0.40 


18 07 30.7 


-19 56 28 


6.5 


(0,-1) 


1110+220 


30 


6 










1900+380 


120 




G12.42+0.50 


18 07 56.4 


-17 56 37 


2.1 


(-10,0) 


210+42 


30 


12 










440+88 


120 




G12.89+0.49 


18 08 56.3 


-17 32 16 


3.5 


(0,0) 


220+44 


30 


8 










340+68 


120 




G12. 21-0.10 


18 09 43.7 


-18 25 09 


13.7 


(-5,5) 


230+46 


30 


2 


G13.87+0.28 


18 11 41.5 


-16 46 34 


4.4 


(10,-8) 


190+38 


30 


1 










430+86 


120 




W33A'* 


18 11 44.0 


-17 53 09 


4.0 


(3,4) 


350+70 


30 


5 










960+190 


120 




G14.33-0.64 


18 16 00.8 


-16 49 06 


2.6 


(0,-4) 


440+88 


30 


8 










830+170 


120 




GL2136^ 


18 19 36.6 


-13 31 40 


2.0 


(8,-9) 


240+48 


30 


9 










520+100 


120 




G19.61-0.23 


18 24 50.1 


-11 58 22 


4.0 


(2,-6) 


460+92 


30 


6 










500+200 


120 
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Table 1 — Continued 



Source 


a (1950.0) 

m s ^ 


S (1950.0) 

(C , 


D 
(kpc) 


Centroid'^ 
("/') 


350/zm 
S." (Jy) 


(") 


D 

Ref. 


G23.95+0.16'' 


18 31 40.8 


-07 57 17 


5.8 


(7,0) 


100±20 


30 


6 












320±64 


120 




G24.49-0.04 


18 33 22.8 


-07 33 54 


3.5 


(7,3) 


190±37 


30 


6 










190±37 


120 




W43S 


18 43 26.7 


-02 42 40 


8.5 


(2,4) 


360±72 


30 


6 










440±88 


120 




G31.41+0.31 


18 44 59.5 


-01 16 07 


7.9 


(-1,0) 


460±92 


30 


1 










740±150 


120 




G40.50+2.54 


18 53 45.6 


07 49 16 


2.1 


(2,1) 


240±48 


30 
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120 






1 O C O C 1 A 
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AO 1 £i on 

Uz lb zy 
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AC riT 1 n 1 Q 
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iU 40 4z 


y. i 
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1 QO 1 9/? 


oU 


/? 



^ A o CI in no 




1 O ACi A A 

\.6 4y 44 


1 1 o 

ii.O 
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(0,-2) 


iUU±zU 


on 













onn 1 /in 
zUUd=4U 


1 on 




tit Tr" 1 TV itH 

W51M° 


19 21 26.2 


-1/1 C\ A O 

14 24 36 


7.0 


(-29,-18) 


280±56 


30 


7 










3690±740 


120 




S87° 


19 44 14.0 


24 28 10 


2.3 


(0,-3) 


140±28 


30 


6 










310±62 


120 




S88B 


19 44 42.0 


25 05 30 


2.0 


(20,-10) 


150±30 


30 


6 










540±110 


120 




ON 1 


20 08 09.9 


31 22 42 


6.0 


(0,-10) 


320±64 


30 


6 










c n 1 i on 

650±130 


120 




ON 2S 


20 19 48.9 


37 15 52 


5.5 


(5,-2) 


200±40 


30 


6 










510±100 


120 




S106^ 


20 25 32.8 


37 12 54 


4.1 


(0,-5) 


110±22 


30 


7 










400±80 


120 




GL2591* 


20 27 35.5 


40 01 13 


1.0 


(6,-2) 


230±46 


30 


3 










440±88 


120 




G97.53+3.19 


21 30 37.0 


55 40 36 


8.5 


(0,-12) 


50±10 


30 


5 












90±18 


120 




BFS 11-B 


21 41 57.6 


65 53 17 


2.0 


(1,-1) 


40±8 


30 


6 










90±18 


120 




S140'' 


22 17 41.1 


63 03 42 


0.90 


(-10,-4) 


350±70 


30 


12 










1210±240 


120 




CEP A 


22 54 19.2 


61 45 44 


0.73 


(-1,-5) 


430±86 


30 


11 










1500±300 


120 




S158 


23 11 36.1 


61 10 30 


2.8 


(-3,-4) 


250±50 


30 


6 










700±140 


120 




NGC 7538-1'^ 


23 11 36.7 


61 11 51 


2.8 


(-1,1) 


150±30 


30 


6 










1240±250 


120 




NGC 7538-9=* 


23 11 52.8 


61 10 59 


2.8 


(3,0) 


130±26 


30 


6 










330±66 


120 




S157 


23 13 53.1 


59 45 18 


2.5 


(-9,16) 


70±14 


30 


6 










280±56 


120 





'^Previously studied by van der Tak et al. 2000. ''Double or multiple peaks in 350 iim map. ^350 fim 
peak offset from the water maser position (0,0). The 350 /im flux density and uncertainty measured in an 
aperture of diameter 6ap- The uncertainty is 20%. For sources with no reported 120" aperture flux density, 
the source was cither very weak (IRAS 02395+6244, IRAS 02461+6147, and G137.07-3.00) or the source 
was very centrally peaked and all of the flux is within a 30" aperture (G12.21— 0.10 and G45.07+0.13). 
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Zinchenko 1994 
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Table 2 
Collected Photometry 
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Table 2 — Continued 



Source 



A 


S 


Beam'' 


Ref. 


(um) 


(Jy) 


(") 




12 


1.9±0.3 


300x45 


1 


25 


11.8±0.7 


300x45 


1 


60 


<189 


90x300 


1 


100 


552±77 


180x300 


1 


350 


30±6 


30 




350 


40±8 


120 




12 


470±14 


300x45 


1 


20 
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25 


20 
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25 
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300x45 


1 


27 
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25 
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25 
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25 
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1 
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25 
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16±0.6 


300x45 


1 


25" 


77±3 


300x45 


1 


60" 


10321±34 


90x300 


1 


100" 


1715±189 


180x300 


1 


350" 


130±26 


30 
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30 




350" 


1580±320 


120 
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Source 


A 


Si/ 


Beam'^ 


Ref. 




(^m) 
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Table 2 — Continued 



Source 
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Beam'^ Ref. 




(^m) 


(Jy) 
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G13.87+0.28 
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40° lOOOilOO 

42° 1300±130 

59° 2350±235 

60 2206±530 

73° 3400±340 

77° 4100±410 

101° 4050±405 

100 6183±1422 

135° 3900±390 

135 4000±400 
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Table 2 — Continued 



Source 


A 




Beam'^ Ref. 




(^m) 


(Jy) 


(") 



G31.41+0.31 



G40.50+2.54 



G35.58-0.03 



G45.07+0.13 



G48.61+0.02 



W51M 



S87 



25" 


1697±136 


300x45 


1 


60" 


7501±525 


90x300 


1 


100" 


11669±3151 


180x300 


1 


350" 


360±72 


30 




350" 


440±88 


120 




1300" 


8 


90 


3 


1300 


20±1.9 


12 


32 


12 


4±0.4 


300x45 


1 


25" 


52±5.2 


300x45 


1 


60" 


1093±197 


90x300 


1 


100" 


2815±394 


180x300 


1 


350 


997±200 


11 


1 


350" 


460±92 


30 
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740±150 


120 




450" 


84±17 


9 


18 


850" 


27±1.4 


15 
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1300 


14 


90 


3 
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4.9±1.0 


22 
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2000" 


2.9±0.6 


34 


18 


12" 


31.8±4.5 


300x45 
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242±24 


300x45 


1 


60" 


2351±423 


90x300 
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100" 


4218±840 


180x300 


1 


350" 


240±48 


30 




350" 


600±120 
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450" 


215±54 


18 


24 


850" 


33±7 


18 


24 


12 


<6.0 


300x45 


1 


25 


77±3.9 


300x45 


1 


60 


1507±196 


90x300 


1 


100 


2594±571 


180x300 


1 


350 


110±22 


30 




350 


120±24 


120 




12 


58±3.5 


300x45 


1 


25 


494±30 


300x45 


1 


60 


<5913 


90x300 


1 


100 


<7497 


180x300 


1 


350 


180±36 


30 




12" 


<25 


300x45 


1 


25" 


175±10 


300x45 


1 


60" 


3195±320 


90x300 


1 


100" 


5227±523 


180x300 


1 


350" 


100±20 


30 




350" 


200±40 


120 




12 


424±42 


300x45 


1 


25 


4344±430 


300x45 


1 


100 


< 26760 


180x300 


1 


350 


280±56 


30 




350 


3690±740 


120 




1300 


27±2.7 


25 


37 


12 


47±3.8 


300x45 


1 


25 


425±26 


300x45 


1 
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Table 2 — Continued 



Source 



A 


S ^ 


Beam'^ 


Ref. 


(um) 


(Jy) 


(") 




60 


3446±310 


90x300 


1 


100 


5158±464 


180x300 


1 


350 


140±28 


30 




350 


310±62 


120 




450 


110±28 


18 


24 


850 


17±3 


18 


24 


1300 


3.6 


90 


3 


10.2 


0.04±0.01 


17 


34 


11.1 


0.12±0.024 


17 


34 


12" 


93±6.5 


300x45 


1 


12.6" 


22±4.4 


17 


34 


17" 


40±8.0 


17 


34 


25" 


1185±71.1 


300x45 


1 


60" 


8686±1129 


90x300 
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100" 


13214±1057 


180x300 


1 


350" 
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3.4 


90 


3 
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300x45 
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90x300 


1 
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3119±312 
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320±64 
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400±80 
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12 
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300x45 


1 


12.6 
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6.6 


28 


19.5" 


630 


6.6 


28 


23" 


920 


6.6 


28 


25 


1112±67 


300x45 


1 
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4600±920 


49 


28 


60 


5314±425 


90x300 


1 


95" 


5800±300 


49 


28 


100 


5721±858 


180x300 


1 


110" 


5500±1100 


49 


28 


160" 


34006±80 


49 


28 


350" 


230±46 


30 





S88B 



ON 1 



ON 2S 



S106 



GL2591 
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Table 2 — Continued 



Source 


A 


Si/ 


Beam^ 


Ref. 




(^m) 


(Jy) 
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oonnj-i Kr\ 
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on 

oU 
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10600ib2650 


20 


11 




60" 


17000zb3400 


90x300 


1 




85 


46500±12900 


270 


27 




100" 


230004ib600 


180x300 


1 




100" 


20200±5050 


30 


11 




125" 


33100±9170 


50 


13 




150 


23400zb6800 


270 


27 




350" 


430zb86 


30 






350" 


1500zb300 


120 
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Table 2 — Continued 



Source 


A 




Beam'' 


Ref. 




(yum) 


(Jy) 


(") 






400° 


2570±741 


50 


13 




450° 


737±140 


20 


33 




550 


27900±6300 


50 


13 




800° 


86±10 


20 


33 




865 


84±18 


8 


16 




1300° 


26±8 


40 


16 




1300 


27±1.6 


30 


45 


S158 


12 


243±7.3 


300x45 


1 




25 


1780±71 


300x45 


1 




60 


7073±495 


90x300 


1 




100 


14138±1131 


180x300 


1 




350 


249±50 


30 


5 




350 


250±50 


30 






350 


700±140 


120 






1300 


15 


30 


5 


NGC 7538-1 


12.5° 


149±21 


7.5 


1 




20° 


250±50 


6 


46 




25° 


640±130 


6 


46 




30° 


2300±700 


40 


46 




50° 


6700±2010 


40 


46 




100° 


11000±3300 


55 


46 




350° 


150±30 


30 






350° 


1240±250 


120 






1000° 


30±9 


55 


46 


NGC 7538-9 


12.5° 


74±13 


9 


46 




20° 


124±30 


6 


46 




25° 


260±50 


6 


46 




30° 


500±150 


40 


46 




50° 


1300±390 


40 


46 




100° 


2700±810 


55 


46 




350° 


130±26 


30 






350° 


330±66 


120 






1000° 


51±5 


55 


46 


S157 


12° 


29±3 


300x45 


1 




25° 


233±12 


300x45 


1 




60° 


1759±123 


90x300 


1 




100° 


264±303 


180x300 


1 




350° 


70±14 


30 






350 


280±56 


120 






850° 


5.9±1.2 


18 


24 



*Flux value used in model. ^Observed flux density and uncertainty. 350 /um flux densities with Beam = 
30" or 120" are from this work. '^FWHM beam size for observed flux density. 

References. — 1. IRAS PSC, 2. Cheung 1980, 3. Chini 1986a, 4. Chini 1986b, 5. Chini 1986c, 6. Chini 
1991, 7. De Buizer 2000, 8. Dent 1988, 9. Dent 1998, 10. Dyck 1977, 11. Ellis 1990, 12. Evans 1979, 
13. Evans 1981a, 14. Evans 1981b, 15. Fazio 1978, 16. Gordon 1990, 17. Guertlcr 1991, 18. Hatchcll 
2000, 19. Hoare 1991, 20. Hefner 1996, 21. Hunter 2000, 22. Jaffc 1984, 23. Jaffa 2001 (priv. comm.), 24. 
Jenness 1995, 25. Jenness 1996, 26. Kastncr 94, 27. Koppcnaal 1979, 28. Lada 1984, 29. Lebofsky 1976, 30. 
McCutcheon 1995, 31. Mooherjea 1999, 32. Moony 1995, 33. Moriarity-Schieven 1991, 34. Pipher 1977, 35. 
Richardson 1989, 36. Sandell 2001, 37. Schoerb 1987, 38. Schwartz 1977, 39. Simon 1985, 40. Soifer 1975, 
41. Sticr 1984, 42. Thronson 1979, 43. Thronson 1980, 44. Tothill 1999, 45. Walker 1990, 46. Werner 1979, 
47. Wright 1977, 48. Zhou 1994, 49. Gezari 1998, 50. Mookerjea 2000 
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Table 3 

Modeled and Observed Source Properties 



Source 


pa 


, 

(10^ cm-3) 


^hof 

(10^ Lq) 


(K) 




IRAS 00338+6312 


1.25 


0.4 


0.11 


46 


1.6 


G123.07-6.31 


1.75 


3.0 


0.60 


51 


1.1 


W3 






2.4 


77 


1.0 


W3(0H) 


1.50 


1.5 


9.5 


65 


1.6 


IRAS 02395+6244 






5.6 


107 


0.8 


IRAS 02461+6147 






1.9 


99 


1.1 


GL490 






0.24 


112 


1.5 


Ori-IRC2 






9.9 


73 


1.5 


S231 


1.50 


1.0 


1.3 


63 


1.5 


S235 






0.95 


125 


2.2 


S241 






1.3 


62 


1.0 


MonR2 






3.7 


98 




S252A 


1.75 


1.0 


0.64 


68 


1.0 


RCW142 


2.25 


6.5 


5.7 


64 


1.0 


W28A2 


2.25 


25 


20 


60 


0.9 


M8E 


1.75 


1.2 


1.7 


93 


1.6 


G9. 62+0. 10 


2.00 


17 


35 


68 


1.1 


G8.67-0.36 


2.00 


15 


13 


67 


1.0 


GlO.60-0.40 


2.50 


120 


92 


52 


1.0 


G12.42+0.50 


2.00 


3.2 


1.9 


70 


1.1 


G12.89+0.49 


2.00 


5.7 


3.9 


56 


1.0 


G12.21-0.10 






55 


57 


0.5 


G13.87+0.28 


1.75 


2.0 


13 


93 


1.6 


W33A 


1.50 


2.0 


10 


65 


1.5 


G14. 33-0.64 


2.00 


1.8 


10 


55 


1.0 


GL2136 


1.75 


6.3 


1.0 


173 


1.1 


G19.61-0.23 






18 


74 


1.1 


G23.95+0.16 


1.50 


0.5 


19 


91 


1.5 


G24.49-0.04 


2.25 


8.5 


4.9 


65 


1.1 


W43S 


2.50 


60 


160 


93 


0.8 


G31.41+0.31 


2.25 


80 


23 


55 


0.8 


G40.50+2.54 


1.50 


0.8 


2.8 


75 


1.6 


G35. 58-0.03 






4.2 


66 


1.3 


G45.07+0.13 






120 


79 


1.3 


G48. 61+0.02 






100 


70 


0.5 


W51M 






380 


88 


1.0 


S87 






4.8 


80 


1.3 


S88B 


1.25 


0.1 


9.0 


75 


2.8 


ONI 


1.75 


7.0 


15 


57 


1.4 


0N2S 


1.75 


2.0 


37 


80 


1.5 


S106 






50 


94 


0.5 


GL2591 


2.00 


0.9 


2.0 


86 


1.1 


BFS 11-B 






0.70 


64 




S140 


1.25 


0.3 


1.9 


91 


2.6 


Cep A 


1.50 


0.5 


2.2 


84 


1.8 


S158 






19 


91 


1.0 


NGC 7538-1 


1.50 


1.2 


15 


74 


1.8 


NGC 7538-9 


1.75 


1.2 


3.6 


81 


1.5 


S157 


0.75 


0.01 


2.9 


86 


3.6 



^The best fit model power law density distribution exponent. ^The model density at 1000 AU. 

•^The bolometric luminosity calculated from the SED. ''The temperature of a blackbody with the 
same mean frequency of the SED. "^The ratio of the deconvolved source size to the beam size. Sources 
with Qdncl^mb < 0.8 were not modeled. 
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Table 4 

Masses axd Sizes 



Source 


(io^'^'au) 




(10^ AU) 


(Mr) 


(Mr) 


Lhol 1 Migo 

(L„/ Mo) 


IRAS 00338+6312 


1.0 


10 


12 


810 


50 


24 


G123.07-6.31 


1.9 


80 


9.4 


570 


220 


27 


W3 


1.6 








180 


140 


W3(0H) 


2.9 


130 


13 


1210 


1010 


94 


IRAS 02395+6244 


4.6 


- 


- 


- 


210 


260 


IRAS 02461+6147 


3.8 


- 


- 


- 


70 


290 


G137.G7-3.00 


1.9 


- 


- 


- 


80 


- 


GL490 


1.0 


- 


- 


- 


20 


100 


Ori-IRC2 


0.5 


- 


- 


- 


240 


410 


S231 


2.7 


70 


9.9 


530 


430 


30 


S235 


2.7 


- 


- 


- 


100 


100 


S241 


3.4 


- 


- 


- 


130 


100 


MonR2 


- 


- 


- 




180 


210 


S252A 


1.2 


10 


5.0 


90 


110 


57 


RCW142 


1.4 


50 


4.8 


130 


410 


140 


W28A2 


1.4 


200 


9.0 


790 


1660 


120 


M8E 


2.2 


40 


5.6 


120 


190 


88 


G9.62+0.10 


4.8 


670 


12 


1660 


5780 


61 


G8.67-0.36 


3.2 


390 


12 


1470 


3650 


36 


GlO.60-0.40 


4.6 


1330 


10 


1960 


12500 


74 


G12.42+0.50 


1.8 


50 


5.4 


140 


300 


64 


G12.89+0.49 


2.5 


120 


7.5 


350 


640 


61 


G12.21-0.10 


5.3 


- 


- 


- 


6770 


81 


G13.87+0.28 


5.3 


190 


7.7 


300 


1290 


100 


W33A 


4.4 


320 


15 


1990 


2390 


42 


G14.33-0.64 


1.9 


150 


9.8 


800 


880 


20 


GL2136 


1.7 


20 


5.2 


90 


320 


200 


G19.61-0.23 


3.4 


- 


- 


- 


1290 


140 


G23.95+0.16 


6.4 


130 


5.9 


120 


1690 


110 


G24.49-0.G4 


3.0 


120 


5.6 


190 


360 


140 


W43S 


4.8 


680 


8.3 


890 


4940 


320 


G31.41+0.31 


4.5 


1500 


15 


3830 


7150 


33 


G40.50+2.54 


2.5 


60 


8.5 


350 


410 


68 


G35. 58-0.03 


3.4 


- 


- 


- 


230 


190 


G45.07+0.13 


9.5 


- 


- 


- 


2620 


460 


G48.61+0.02 


4.6 


- 


- 


- 


4350 


230 


W51M 


5.0 


- 


- 


- 


28100 


140 


S87 


2.2 


- 


- 


- 


260 


190 


S88B 


4.4 


50 


4.6 


50 


260 


270 


ONI 


6.3 


810 


15 


2390 


3650 


37 


0N2S 


6.0 


220 


7.2 


270 


2410 


150 


S106 


0.2 


- 


- 


- 


1030 


490 


GL2591 


0.8 


10 


2.9 


20 


68 


290 


t .Oo-t-o.iy 










iUoU 




BFS 11-B 










60 


130 


S140 


1.8 


20 


8.0 


250 


150 


120 


Cep A 


1.0 


10 


11 


280 


120 


180 


S158 


2.0 








860 


220 


NGC 7538-1 


3.7 


150 


11 


750 


1510 


100 


NGC 7538-9 


3.1 


60 


5.6 


120 


400 


180 


S157 


6.8 


50 


3.9 


10 


270 


110 



"^The radius of the deconvolved FWHM size, ^deo of the 350 /zm emission. ^The integrated mass within 
rdec- "The density equals 10"* cm~^ at this radius. ''The integrated mass within r„. ''Calculated with the 
isothermal temperature from Eq. (6) for the modeled sample and (Tj^o) = 29 K for the sources that were 
not modeled. 



